Antiferromagnetic Spatial Ordering in a Quenched One-dimensional Spinor Gas 
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We have experimentally observed the emergence of spontaneous antiferromagnetic spatial order in 
a sodium spinor Bose-Einstein condensate that was quenched through a magnetic phase transition. 
For negative values of the quadratic Zeeman shift, a gas initially prepared in the F = l,mj? = 
state collapsed into a dynamically evolving superposition of all 3 spin projections, rar = 0, ±1. The 
quench gave rise to rich, nonequilibrium behavior where both nematic and magnetic spin waves were 
generated. We characterized the spatiotemporal evolution through two particle correlations between 
atoms in each pair of spin states. These revealed dramatic differences between the dynamics of the 
spin correlations and those of the spin populations. 

PACS numbers: 03.75.Mn,67.85.De,67.85.Fg,67.85.Hj 



Dynamics far from equilibrium occur ubiquitously in 
nature. A prime example can be found close to a symme- 
try breaking phase transition. A rapid passage through 
the critical point can quickly nucleate spatial inhomo- 
geneities that subsequently evolve slowly with time. Such 
behavior typifies physical systems of vastly different mi- 
croscopic origin, for example, rapid cooling of the early 
universe following inflation [l|-|3|, the chiral symmetry 
breaking phase transition in heavy-ion collisions and 
sudden cooling of solid-state magnets below the Curie 
temperature |6j. The nature of the pattern formation is 
of fundamental interest to uncover its universal features 
0, El- For example, the temporal evolution of correla- 
tion functions can yield critical exponents of the phase 
transition [§] . Understanding how entanglement and cor- 
relations flow in a many-body system is also relevant for 
enabling quantum computation [Toj . 

The spin degrees of freedom of a Bose gas offer exciting 
possibilities for the exploration of non-equilibrium phe- 
nomena In this work we have realized a quantum 
quench in the laboratory by crossing through a magnetic 
phase transition in an antiferromagnetic, F = 1 sodium 
spinor Bose-Einstein condensate (BEC), as shown in Fig- 
ure [TJ We have captured the stochastic dynamics of this 
quench through spatially resolved measurements of two 
particle spin correlations. The spin-dependent mean-field 
Hamiltonian, 

H sp = c 2 n(F) 2 + q(F 2 z ) 

contains a competition between intrinsic magnetic in- 
teractions and the coupling to an external magnetic 
field. The latter is parameterized by the quadratic Zee- 
man shift q, which we rapidly switched from positive 
to negative values to initiate the quench. Similar cor- 
relation measurements have been performed recently for 
quenched superfluids without spin [3, EH ■ In the above 
formula, n is the mean particle density and ci is the 
strength of spin-dependent interactions, which can be ei- 
ther positive (for antiferromagnets) or negative (for ferro- 
magnets). F and F z are the vector spin 1 operator and its 



z projection, respectively. A variety of quantum phases 
can be realized depending on the value of q 11] . In F = 1 
87 Rb, where c<i < 0, altering q suddenly through a phase 
boundary resulted in spontaneous transverse magnetiza- 
tion patterns (l6| . 
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FIG. 1: (Color online). Quantum quench dynamics, (a) 
Instantaneous change in quadratic Zeeman energy q brings 
the gas through a phase transition between polar (ttif = 0) 
and antiferromagnetic (m? = ±1) ground states at t = 0. 
Elongated BEC (yellow ellipse) subsequently develops one- 
dimensional spin domains, (b) Time-of-flight Stern-Gerlach 
images show the spontaneous axial spin structures at a hold 
time of 48 ms (upper image) while by 2.5 s (lower image) the 
system has relaxed to an antiferromagnetic configuration free 
of domains. Residual fringes in the latter image are experi- 
mental artifacts due to the probe light. 

The case C2 > has been much less studied. Here it 
is nematic, rather than ma gne tic, order that is the prin- 
cipal degree of freedom 17H20] , which profoundly affects 
the nature of the spatial patterns and their dynamics 
21]. Early work on antiferromagnetic gases principally 
explored instabilities in two spin components, where ne- 
matic order played no role [221 ] . In our experiment, we 
observe spatial dynamics after tuning q through a phase 
transition at q = using microwave magnetic fields, a 
transition first observed in our earlier work (23| . Related 
work on F — 2 87 Rb BEC, for which C2 > 0, triggered 
pattern formation at a specific wavelength by preparing 
an unstable initial state [24j . Our work explores comple- 
mentary dynamics, where passage through a phase tran- 
sition generates fluctuations in the nematic order on all 
length scales rather than just a specific wavevector as in 
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24J . Thus the path to equilibrium involves the collective 
dynamics of many modes, which is our focus. 

We used time-of-flight Stern-Gerlach (TOF-SG) obser- 
vation to record the one-dimensional spatiotemporal pat- 
tern formation in all 3 spin components, mp = 0, ±1, 
with high resolution. This method works only for ex- 
treme aspect ratio traps (25[. Importantly, we also re- 
solved the spin density fluctuations, which revealed the 
crucial role played by the spin-mixing interaction be- 
tween atom pairs, 2(mf = 0) <-> (mp = 1) + (mp = — 1). 
This interaction leads to entanglement between spins 26- 

, and to density correlations since the fluctuations are 
simply related by Auq = — 2An_i = — 2An + i. In our ex- 
periment we show these to initially develop locally, and 
then to evolve globally. We have spatially resolved the 
full spin-density correlation matrix of a quantum quench, 
to our knowledge for the first time. In thermal equi- 
librium these density correlations directly determine the 
compressibility matrix (29L l30| . 

Optically confined Bose-Einstein condensates in the 
mp = state were prepared in a static magnetic field 
of B x = 100 mG in a manner described in our earlier 
work [23j . Transverse magnetic fields and field gradients 
were compensated to within 5 mG and 0.6 mG/cm, re- 
spectively. The peak density uq = 5 x 10 14 cm~ 3 and 
axial Thomas-Fermi radius R x — 340 /im were measured 
to an accuracy of 5%, from which we determined the 
spin-dependent interaction energy c 2 no = h x 120 Hz. 
The axial and radial trapping frequencies were 7 and 470 
Hz, respectively, accurate to 10%. The radial Thomas- 
Fermi radius, R,± = 5 /im, was small enough such that 
only axial spin domains could form pjjj . The measured 
temperature was 400 nK, close to the chemical potential 
of 360 nK. 

We rapidly switched q from qi = h x +2.8 Hz > to 
a final value qf = h x —4.2 Hz < at t = 0. A quantum 
phase transition divides the polar (mp — 0) ground state 
from an "antiferromagnetic" one (mp = ±1 superposi- 
tion) at q = 0, as shown in Figure [TJ The dynamical 
evolution from one to the other state is governed by the 
spin-dependent interaction term oc c 2 . q was switched by 
turning on a far-off resonance microwave "dressing" field 
tuned to near the F = 1 — > F = 2 ground state hyper- 
fine resonance. The value of qf was controlled through 
the AC Stark shift by adjusting the microwave power 

& 

Since qi, \qf\ -C c 2 ng, the gas remained very close to 
the phase transition point at all times. Nonetheless the 
change triggered a rapid instability in the mp = con- 
densate. Through spin exchange collisional interactions 
we observed that the mp = ±1 population fraction, f pm , 
increased rapidly near t = 20 ms (see Figure [4jo) . This 
was accompanied by oscillations and a slower period of 
relaxation over > 100 ms toward an apparent equilib- 
rium. We observed a 30% atom loss over 100 ms due to 
off-resonant microwave excitation. 



After holding the gas for a time t, the initial spin den- 
sity fluctuations were amplified, and ultimately, macro- 
scopic, one-dimensional spatial domains were observed to 
have formed (see FigureQ})). We then diagnosed the axial 
spin distribution using a TOF-SG sequence [32j consist- 
ing of 2 ms TOF, 4 ms pulsed SG field and 24 ms addi- 
tional TOF to separate the 3 spin states, followed by a 50 
fis absorption imaging pulse on the F = 2 — > F' = 3 tran- 
sition of the D2 line. Due to the extreme 70:1 aspect ra- 
tio the TOF axial distribution remained very close to its 
in-situ value. Each image was Fourier filtered to remove 
spurious interference fringes with wavelength between 60- 
62 /im created by the probe light, and summed over the 
radial direction to obtain one-dimensional atom density 
distributions. We took the center of these distributions 
(x = 0) to be the maximum of the average density profile 
for 30 runs at each hold time t, after smoothing by a 200 
/im moving average filter. 
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FIG. 2: (Color online.) Sample profiles no(x) and n_i(x) 
measured on a single shot at t = 68 ms. 

Figure [5] shows measured spin density profiles, uq(x) 
and n-i(x), for a single shot of the experiment at a hold 
time of 68 ms. Domains as small as 13 /im half-width at 
half-maximum can be seen. This is 4% of the Thomas- 
Fermi radius, and is comparable to our imaging resolu- 
tion, d w 10 /im [39]. Domains smaller than d would have 
expanded to a size > d, and thus cannot be distinguished 
by the time-of-flight method from larger scale features. 
The energetics of the quench, however, suggest that do- 
mains no smaller than h/ (2y / Mqf) = 5/im are likely to 
appear. 

While single images contained multiple spin domains, 
the average over many runs showed no structure, indi- 
cating that the formation and dynamics were stochas- 
tic in nature. Therefore, we quantified the data using 
two-particle correlations, specifically computing the two- 
point covariance function (TPCF) for the density fluctu- 
ations An of a spin pair i, j, 

S ilj (x 1 ,X2,t) = (An i (xi,t)Ah j (x 2 ,t)) 

where (...) refer to ensemble averages of fluctuations 
about the sample mean, taken over 30 runs of the ex- 
periment. Averages over a smaller sample size yielded 
similar results. We further analyzed these data in terms 
of both the local correlation function (x\ = x 2 ) as well as 
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non-local effects (x± ^ x^). We discuss both these results 
in turn below. 



as 

cu 1 — 1 



as as 

(J cu 

O D_ 



-E- R [-1,-1] 
-■- R [+1,-1] 
-1/2 R [0,-1] 








20 








•elati 
[art 




* 








Local Ct 
Peak F 







40 



60 



80 100 



Hold time after quench [ms] 



--B-- R [-1,-1] 
R [+1,-1] 
-1/2 R [0,-1] 



100 200 300 400 

Hold time after quench [ms] 



500 



FIG. 3: (Color online.) Correlation of atom pairs after the 
quench at t = 0. (a) Measured local density correlations 
— J dx(Afii(x,t)Ahj(x,t)} for spin states (i, j) = 
( — 1, —1) and (+1, —1) versus time. Also shown is the (0, —1) 
correlation scaled by a factor — i as described in the text, (b) 
Correlation functions determined from numerical simulations. 
The horizontal dashed line marks zero correlation. 



Figure [3^ shows the measured time evolution 
of the local pair correlation function, Rij{t) — 
J dx{Ani(x,t)Anj(x,t)). This graph is a central result 
of this work. Noise due to total atom number fluctua- 
tions was subtracted from the data as described later in 
the manuscript. For short times, t < 35 ms, all 3 spin 
states exhibit correlations whose absolute value rises from 
zero to a maximum value. This is a result of the local 
spin mixing process discussed earlier, which predicts that 
the ratio Rx -i/R-i _i = 1.0 and R ,-i/R-i ,-i = -2.0. 
Our measured peak values of the ratios are 0.9 and —2.0, 
respectively. Thus our data indicate that although the 
system is outside the single-mode regime 



11], spin mix- 



ing generates nearly the maximum local correlations that 
are possible. 

What is equally interesting is that for t > 35 ms the 
data all diverge from one another. The diagonal corre- 
lations, appear to reach a steady state of 50 % 
of their peak value. This is the result of saturation of 
the parametric gain due to the finite number of rriF = 
atoms, which causes fluctuations in the amplified modes 
to be suppressed By contrast, the off-diagonal data, 
decay sharply. The i?o,-i data has been scaled 
by the factor — i to allow for comparison of all 3 curves, 
and exhibits an intermediate behavior between 
and 

The initial, positive correlation between ±1 atoms 



is indicative of nematic spin wave excitations, which 
are small oscillations in the order parameter N zz = 
n + i(x,t) + n_i(x,t) — 2n Q (x,t). These excitations oc- 
cur naturally for F = 1 spinor condensates with c~2 > 0, 
since mp — regions phase separate from mp = ±1 [32] ]. 
As time evolves, distinct regions in the condensate can 
communicate with one another through the propagation 
of these spin waves. The decay of R+\^-\ can therefore be 
viewed as a destructive intererence between waves origi- 
nating from different points in space. Interference effects 
have also been recently reported in quenched scalar su- 
perfluids, leading to oscillations in the density-density 
correlation function jl5j |. 

Accompanying this decay of J2+i,_i is a growth of 
magnetic spin wave excitations, whose order param- 
eter is M z (x,t) = n + i(x,t) — ri-i(x,t). The local 
magnetization variance is vax[M z (x, t)] = S+1,+1 + 
S-i^—i — 2S'_i,+i, which becomes non-zero for imper- 
fect correlation between the two spins, S—i t +i(x,t) < 
+1 (x, t), S—i t -i(x, t), as we observed experimentally. 
Our data in Figure [3K show that R + i_i reaches zero at 
~ 85 ms, after which the two spins even become nega- 
tively correlated. 

In order to understand our data theoretically, we have 
performed numerical simulations of the 3 coupled spinor 
Gross-Pitaevskii (GP) equations. Assuming a BEC ini- 
tially at zero temperature, we obtained the initial wave- 
function for the mp = component numerically. Vac- 
uum noise in the =±1 states was simulated as clas- 
sical noise according to the Truncated Wigner approxi- 
mation (TWA) [33rl36(. Thirty separate simulations al- 
lowed us to compute the expectation value of quantum 
mechanical observables as statistical averages over dif- 
ferent random initial conditions. Vacuum modes with 
wavelength less than t; sp i n are not expected to contribute 
to the spin instability. Therefore, we imposed a cutoff 
energy of C2Uq, which resulted in ~700 virtual particles, 
while the condensate contained 5 x 10 6 particles, similar 
to the experimental conditions. The results were rela- 
tively insensitive to this cutoff as expected for an expo- 
nential amplification process. For example, magnifying 
the noise by a factor of 10 changed the time at which the 
±1 fraction = 0.5 by only 20 %. 

Figure shows the numerically obtained correlation 
functions, which shows good qualitative agreement with 
the experimentally observed trends. Similar to the ex- 
perimental data, R+i t -i has dropped below zero in the 
simulation after sa 300 ms. The overall dynamics were 
faster in the experiment by approximately this factor, 
300/85 = 3.5, which we cannot presently account for. 
Possible effects include residual ±1 population in the ini- 
tial state caused by thermal or technical fluctuations that 
would speed up the instability. In addition, atom loss (in 
our case due to off-resonant microwave excitation) can 
cause decoherence (28j . 

The simulations show that even in the absence of ther- 
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mal or technical noise, vacuum fluctuations can induce 
a decay of correlations, a somewhat counterintuitive re- 
sult. For any particular realization of the experiment, 
the vacuum noise of the two bosonic fields, ip + i and V>_i, 
is uncorrelated. Tiny differences in the initial conditions 
become amplified by the spin mixing instability, resulting 
in final wavefunctions tp±i(t — tp) that differ macroscop- 
ically, i.e., by much more than just the vacuum noise. By 
contrast, identical initial conditions for ip±i (an unphys- 
ical situation) resulted in identical time evolution and no 
decay, consistent with the symmetry of the GP equations 
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FIG. 4: (Color online.) Nonlocal correlations and spin domain 
dynamics, (a) Top row shows Tij(S) = J Sij(x,x + S)dx 
for various spin pairs and relative position 8. Bottom row 
shows the same data with the background caused by atom 
number variations subtracted, (b) Dynamics of the fraction of 
atoms in ±1 (green squares, right axis), the RMS momentum 
width determined from the Fourier transform of T_i._i and 
T-1,+1 (red and blue circles, respectively, left axis) and from 
a theoretical prediction (solid red line, left axis). 



oscillations. The height of the peak was determined by 
an additional Gaussian fit, and constitutes the data of 
Figure [3] 

The oscillations signify a dynamically evolving non- 
local order induced by the rapid quench, which we charac- 
terized by computing the Fourier transforms, T_i._i(/c) 
and T+i i _i(fc). These showed an increase in high spatial 
frequency components with time after the quench. Figure 
|4)d shows this trend as a growth in the second moment of 
the Fourier transforms, k rms = (f^° k 2 Tdk/ J °° Td/c) 1 ^ 2 . 
Higher order effects in the SG expansion are responsible 
for the slight magnification of k rms for the (—1, —1) rel- 
ative to (+1, —1). 

While a rapid quench is expected to lead to coarsening 
of the domains on long timescales @, [37j , our data show 
the opposite trend within the time window 0-100 ms. 
This suggests that it might be an intermediate phase of 
evolution that precedes coarsening. In fact, energy con- 
servation requires the kinetic energy to grow with the 
nip = ±1 populations, since they reduce the quadratic 
Zeeman energy. Using h 2 k 2 ms / (2M) from the correlation 
functions as an estimate of the mean kinetic energy, we 
obtain - ^/-2Mq f f ±1 /h 2 = f ±l (tf/ 2 /(47pm). 

This formula is plotted as a solid line Figure HJd using our 
measurements, and shows the same trend as our data. It 
consistently lie above our measurements, most likely due 
to the fact that we do not account for variations in the 
phase of the order parameter nor the interaction energy 
contained in the finite width domain boundaries. 

In summary, we have measured the time evolution of 
spin correlations in a quenched spinor BEC. Our data re- 
veal rich non-equilibrium dynamics where both nematic 
and magnetic spin density waves are spontaneously gen- 
erated. Future work will explore possible universal sig- 
natures in the decay of correlations [ij]. 
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We conclude by quantifying the non-local dynam- 
ics. Figure 2] shows a one-dimensional reduction of the 
covariance map Sij in terms of relative coordinate 5, 
Tij(5,t) — J(Am(x, t)Arij(x + S, t))dx, for spin pairs 
(-1,-1), (+1,-1) and (0,-1), at t =20 ms after the 
quench. A narrow peak, either positive or negative, can 
be seen near 6 = 0. A smooth, positive background lies 
underneath it, primarily caused by atom number varia- 
tions from one shot to the next. The latter lead to cor- 
related density fluctuations throughout the cloud 



41]. 



The wings of the data appeared to fit well to a Gaus- 
sian distribution. When subtracted, the data revealed 
the intrinsic fluctuations, which are plotted underneath, 
showing a sharp, local correlation peak surrounded by 
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